Nonlinear dynamics of river runoff elucidated by horizontal visibility graphs

Horizontal Visibility Graphs (HVGs) are a recently developed method to construct networks from time series. The values of the time series are considered as the nodes of the network and are linked to each other if there is no larger value between them, such as they can “see” each other. The network properties reflect the nonlinear dynamics of the time series. For some classes of stochastic processes and for periodic time series, analytical results can be obtained for network-derived quantities such as the degree distribution, the local clustering coefficient distribution, the mean path length, and others. HVGs have the potential to discern between deterministic-chaotic and correlated-stochastic time series. Here, we investigate the sensitivity of the HVG methodology to properties and pre-processing of real-world data, i.e., time series length, the presence of ties, and deseasonalization, using a set of around 150 runoff time series from managed rivers at daily resolution from Brazil with an average length of 65 years. We show that an application of HVGs on real-world time series requires a careful consideration of data pre-processing steps and analysis methodology before robust results and interpretations can be obtained. For example, one recent analysis of the degree distribution of runoff records reported pronounced sub-exponential “long-tailed” behavior of North American rivers, whereas another study of South American rivers showed hyper-exponential “short-tailed” behavior resembling correlated noise. We demonstrate, using the dataset of Brazilian rivers, that these apparently contradictory results can be reconciled by minor differences in data-preprocessing (here: small differences in subtracting the seasonal cycle). Hence, data-preprocessing that is conventional in hydrology (“deseasonalization”) changes long-term correlations and the overall runoff dynamics substantially, and we present empirical consequences and extensive simulations to investigate these issues from a HVG methodological perspective. After carefully accounting for these methodological aspects, the HVG analysis reveals that the river runoff dataset shows indeed complex behavior that appears to stem from a superposition of short-term correlated noise and “long-tailed behaviour,” i.e., highly connected nodes. Moreover, the construction of a dam along a river tends to increase short-term correlations in runoff series. In summary, the present study illustrates the (often substantial) effects of methodological and data-preprocessing choices for the interpretation of river runoff dynamics in the HVG framework and its general applicability for real-world time series.Fil: Lange, Holger. Norwegian Institute of Bioeconomy Research; NoruegaFil: Sippel, Sebastian. Norwegian Institute of Bioeconomy Research; NoruegaFil: Rosso, Osvaldo Aníbal. Instituto Universidad Escuela de Medicina del Hospital Italiano; Argentina. Universidad de Los Andes; Chile. Universidade Federal de Alagoas; Brasil. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentin

Robust detection and attribution of climate change under interventions

Fingerprints are key tools in climate change detection and attribution (D&A) that are used to determine whether changes in observations are different from internal climate variability (detection), and whether observed changes can be assigned to specific external drivers (attribution). We propose a direct D&A approach based on supervised learning to extract fingerprints that lead to robust predictions under relevant interventions on exogenous variables, i.e., climate drivers other than the target. We employ anchor regression, a distributionally-robust statistical learning method inspired by causal inference that extrapolates well to perturbed data under the interventions considered. The residuals from the prediction achieve either uncorrelatedness or mean independence with the exogenous variables, thus guaranteeing robustness. We define D&A as a unified hypothesis testing framework that relies on the same statistical model but uses different targets and test statistics. In the experiments, we first show that the CO2 forcing can be robustly predicted from temperature spatial patterns under strong interventions on the solar forcing. Second, we illustrate attribution to the greenhouse gases and aerosols while protecting against interventions on the aerosols and CO2 forcing, respectively. Our study shows that incorporating robustness constraints against relevant interventions may significantly benefit detection and attribution of climate change

Drought, Heat, and the Carbon Cycle: a Review

Purpose of the Review Weather and climate extremes substantially affect global- and regional-scale carbon (C) cycling, and thus spatially or temporally extended climatic extreme events jeopardize terrestrial ecosystem carbon sequestration. We illustrate the relevance of drought and/or heat events (“DHE”) for the carbon cycle and highlight underlying concepts and complex impact mechanisms. We review recent results, discuss current research needs and emerging research topics. Recent Findings Our review covers topics critical to understanding, attributing and predicting the effects of DHE on the terrestrial carbon cycle: (1) ecophysiological impact mechanisms and mediating factors, (2) the role of timing, duration and dynamical effects through which DHE impacts on regional-scale carbon cycling are either attenuated or enhanced, and (3) large-scale atmospheric conditions under which DHE are likely to unfold and to affect the terrestrial carbon cycle. Recent research thus shows the need to view these events in a broader spatial and temporal perspective that extends assessments beyond local and concurrent C cycle impacts of DHE. Summary Novel data streams, model (ensemble) simulations, and analyses allow to better understand carbon cycle impacts not only in response to their proximate drivers (drought, heat, etc.) but also attributing them to underlying changes in drivers and large-scale atmospheric conditions. These attribution-type analyses increasingly address and disentangle various sequences or dynamical interactions of events and their impacts, including compensating or amplifying effects on terrestrial carbon cycling.publishedVersio

Insight into hydrochemistry: a multi-catchment comparison using Horizontal Visibility Graphs

Long time series of environmental variables are reflecting the dynamics and functioning of ecosystems. Here, we investigate data from a long-term monitoring site in Germany, the Bramke valley in the Harz mountains, where time series of ion concentrations in stream water are obtained since the 1970ies at four measurement locations from three small adjacent forested catchments. Since for (only) one of the catchments daily runoff rates are also available, we invent a method to generate time series of nutrient output from the catchments. Both concentrations and outputs show a number of remarkable long-term changes, including ones not obviously related to changes in atmospheric deposition, management or properties of the forest stands. For the analysis of the Bramke data, we investigate Horizontal Visibility Graphs (HVGs), a recently developed method to construct networks based on time series. Values (the nodes of the network) of the time series are linked to each other if there is no value higher between them. The network properties, such as the degree and distance distributions, reflect the nonlinear dynamics of the time series. For certain classes of stochastic processes and for periodic time series, analytic results can be obtained for some network properties. HVGs have the potential to discern between deterministic-chaotic and correlated-stochastic time series. We classify the Bramke series according to their stochastic nature, with a focus on inter-catchment comparison on one hand, on different nutrients for one catchment on the other, and conclude on possible reasons for the observed changes and their ecological interpretation

Have precipitation extremes and annual totals been increasing in the world’s dry regions over the last 60 years?

Daily precipitation extremes and annual totals have increased in large parts of the global land area over the past decades. These observations are consistent with theoretical considerations of a warming climate. However, until recently these trends have not been shown to consistently affect dry regions over land. A recent study, published by Donat et al. (2016), now identified significant increases in annual-maximum daily extreme precipitation (Rx1d) and annual precipitation totals (PRCPTOT) in dry regions. Here, we revisit the applied methods and explore the sensitivity of changes in precipitation extremes and annual totals to alternative choices of defining a dry region (i.e. in terms of aridity as opposed to precipitation characteristics alone). We find that (a) statistical artifacts introduced by data pre-processing based on a time-invariant reference period lead to an overestimation of the reported trends by up to 40 %, and that (b) the reported trends of globally aggregated extremes and annual totals are highly sensitive to the definition of a "dry region of the globe". For example, using the same observational dataset, accounting for the statistical artifacts, and based on different aridity-based dryness definitions, we find a reduction in the positive trend of Rx1d from the originally reported +1.6 % decade−1 to +0.2 to +0.9 % decade−1 (period changes for 1981–2010 averages relative to 1951–1980 are reduced to −1.32 to +0.97 % as opposed to +4.85 % in the original study). If we include additional but less homogenized data to cover larger regions, the global trend increases slightly (Rx1d: +0.4 to +1.1 % decade−1), and in this case we can indeed confirm (partly) significant increases in Rx1d. However, these globally aggregated estimates remain uncertain as considerable gaps in long-term observations in the Earth's arid and semi-arid regions remain. In summary, adequate data pre-processing and accounting for uncertainties regarding the definition of dryness are crucial to the quantification of spatially aggregated trends in precipitation extremes in the world's dry regions. In view of the high relevance of the question to many potentially affected stakeholders, we call for a well-reflected choice of specific data processing methods and the inclusion of alternative dryness definitions to guarantee that communicated results related to climate change be robust.Have precipitation extremes and annual totals been increasing in the world’s dry regions over the last 60 years?publishedVersio

Detecting impacts of extreme events with ecological in situ monitoring networks

Extreme hydrometeorological conditions typically impact ecophysiological processes on land. Satellite-based observations of the terrestrial biosphere provide an important reference for detecting and describing the spatiotemporal development of such events. However, in-depth investigations of ecological processes during extreme events require additional in situ observations. The question is whether the density of existing ecological in situ networks is sufficient for analysing the impact of extreme events, and what are expected event detection rates of ecological in situ networks of a given size. To assess these issues, we build a baseline of extreme reductions in the fraction of absorbed photosynthetically active radiation (FAPAR), identified by a new event detection method tailored to identify extremes of regional relevance. We then investigate the event detection success rates of hypothetical networks of varying sizes. Our results show that large extremes can be reliably detected with relatively small networks, but also reveal a linear decay of detection probabilities towards smaller extreme events in log–log space. For instance, networks with  ≈  100 randomly placed sites in Europe yield a  ≥  90 % chance of detecting the eight largest (typically very large) extreme events; but only a  ≥  50 % chance of capturing the 39 largest events. These findings are consistent with probability-theoretic considerations, but the slopes of the decay rates deviate due to temporal autocorrelation and the exact implementation of the extreme event detection algorithm. Using the examples of AmeriFlux and NEON, we then investigate to what degree ecological in situ networks can capture extreme events of a given size. Consistent with our theoretical considerations, we find that today's systematically designed networks (i.e. NEON) reliably detect the largest extremes, but that the extreme event detection rates are not higher than would be achieved by randomly designed networks. Spatio-temporal expansions of ecological in situ monitoring networks should carefully consider the size distribution characteristics of extreme events if the aim is also to monitor the impacts of such events in the terrestrial biosphere

ANALYSIS OF VTVL AND VTHL REUSABLE LAUNCH VEHICLE CONFIGURATIONS

The German Aerospace Centre (DLR) is conducting systematic analyses of reusable space transportation configurations. Two-stage vertical take-off vertical landing (VTVL) and winged, vertical take-off horizontal landing (VTHL) partially reusable launcher configurations are systematically analyzed. The investigated configurations consider reusable first stages that either perform a return to launch site or land down range of the launch site. The propellant combinations analyzed include LOX/LH2, LOX/LCH4 and LOX/RP-1. Staged combustion and gas generator cycle engines are taken into account. The same type of engines with different expansion ratios are used on the reusable first stages and the expendable upper stages. Major analysis objectives are the comparison of various reusable launch vehicle configurations under similar design assumptions as well as the identification of their critical aspects, benefits and drawbacks

Reverse engineering model structures for soil and ecosystem respiration: the potential of gene expression programming

Accurate model representation of land-atmosphere carbon fluxes is essential for climate projections. However, the exact responses of carbon cycle processes to climatic drivers often remain uncertain. Presently, knowledge derived from experiments, complemented with a steadily evolving body of mechanistic theory provides the main basis for developing such models. The strongly increasing availability of measurements may facilitate new ways of identifying suitable model structures using machine learning. Here, we explore the potential of gene expression programming (GEP) to derive relevant model formulations based solely on the signals present in data by automatically applying various mathematical transformations to potential predictors and repeatedly evolving the resulting model structures. In contrast to most other machine learning regression techniques, the GEP approach generates "readable" models that allow for prediction and possibly for interpretation. Our study is based on two cases: artificially generated data and real observations. Simulations based on artificial data show that GEP is successful in identifying prescribed functions with the prediction capacity of the models comparable to four state-of-the-art machine learning methods (Random Forests, Support Vector Machines, Artificial Neural Networks, and Kernel Ridge Regressions). Based on real observations we explore the responses of the different components of terrestrial respiration at an oak forest in south-east England. We find that the GEP retrieved models are often better in prediction than some established respiration models. Based on their structures, we find previously unconsidered exponential dependencies of respiration on seasonal ecosystem carbon assimilation and water dynamics. We noticed that the GEP models are only partly portable across respiration components; the identification of a "general" terrestrial respiration model possibly prevented by equifinality issues. Overall, GEP is a promising tool for uncovering new model structures for terrestrial ecology in the data rich era, complementing more traditional modelling approaches

Contrasting biosphere responses to hydrometeorological extremes: revisiting the 2010 western Russian heatwave

Combined droughts and heatwaves are among those compound extreme events that induce severe impacts on the terrestrial biosphere and human health. A record breaking hot and dry compound event hit western Russia in summer 2010 (Russian heatwave, RHW). Events of this kind are relevant from a hydrometeorological perspective, but are also interesting from a biospheric point of view because of their impacts on ecosystems, e.g., reductions in the terrestrial carbon storage. Integrating both perspectives might facilitate our knowledge about the RHW. We revisit the RHW from both a biospheric and a hydrometeorological perspective. We apply a recently developed multivariate anomaly detection approach to a set of hydrometeorological variables, and then to multiple biospheric variables relevant to describe the RHW. One main finding is that the extreme event identified in the hydrometeorological variables leads to multidirectional responses in biospheric variables, e.g., positive and negative anomalies in gross primary production (GPP). In particular, the region of reduced summer ecosystem production does not match the area identified as extreme in the hydrometeorological variables. The reason is that forest-dominated ecosystems in the higher latitudes respond with unusually high productivity to the RHW. Furthermore, the RHW was preceded by an anomalously warm spring, which leads annually integrated to a partial compensation of 54% (36% in the preceding spring, 18% in summer) of the reduced GPP in southern agriculturally dominated ecosystems. Our results show that an ecosystem-specific and multivariate perspective on extreme events can reveal multiple facets of extreme events by simultaneously integrating several data streams irrespective of impact direction and the variables' domain. Our study exemplifies the need for robust multivariate analytic approaches to detect extreme events in both hydrometeorological conditions and associated biosphere responses to fully characterize the effects of extremes, including possible compensatory effects in space and time